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The effective equation of motion is derived for a scalar field interacting with other fields in a 
Friedman-Robertson-Walker background space-time. The dissipative behavior reflected in this ef- 
fective evolution equation is studied both in simplified approximations as well as numerically. The 
relevance of our results to inflation are considered both in terms of the evolution of the inflaton field 
as well as its fluctuation spectrum. A brief examination also is made of supersymmetric models that 
yield dissipative effects during inflation. 

PACS numbers: 98.80.Cq 

in 

o 

Q ■ Version published: Physical Review D71, 023513 (2005) 

(N 

§ ■ I. INTRODUCTION 

Hi ! 

It is by now well established that inflation models in general have dissipative effects during the inflationary period 
HI L4 Lj. [H LJ These effects have two major consequences for inflationary dynamics. First they result in radiation 
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production during inflation, which in turn influences the fluctuations that seed large scale structure. Second they lead 
to temporally nonlocal terms in the effective evolution equation of the inflaton, which can significantly influence the 
nature and history of the inflation period. The realization of these dissipative effects in inflation models has resulted 

^C) . in the division of inflation into two dynamical possibilities referred to as cold and warm inflation. Cold inflation is 
simply the original picture of inflation envisioned in the earliest works ja, 0, M, ■ In this picture the effects of dissipation 

^»0 t are negligible. The fluctuations created during inflation arc effectively zero-point ground state fluctuations and the 
evolution of the inflaton field is governed by a ground state evolution equation. In contrast, in the warm inflation 
picture |9(, inflationary expansion and radiation production occur concurrently. In this picture, the fluctuations 
created during inflation emerge from some excited statistical state and the evolution of the inflaton has dissipative 

r^ terms arising from the interaction of the inflaton with other fields. 

The dividing point between warm and cold inflation is roughly at p r ~ H, where p r is the radiation energy density 
Qh! present during inflation and H is the Hubble parameter. Thus p T > H is the warm inflation regime and p r < H 
is the cold inflation regime. This criteria is independent of thcrmalization, but if such were to occur, one sees this 
criteria basically amounts to the warm inflation regime corresponding to when T > H. This is easy to understand 
since the typical inflaton mass during inflation is m^ w H and so when T > H, thermal fluctuations of the inflaton 
field will become important. This criteria for entering the warm inflation regime turns out to require the dissipation 
of a very tiny fraction of the inflaton vacuum energy during inflation. For example, for inflation with vacuum (i.e. 
potential) energy at the GUT scale ~ 10 15_16 GeV, in order to produce radiation at the scale of the Hubble parameter, 
which is ps 10 10_11 GeV, it just requires dissipating one part in 10 20 of this vacuum energy density into radiation. 
Thus e nerg etically not a very significant amount of radiation production is required to move into the warm inflation 
regime [H|. In fact the levels are so small, and their eventual effects on density perturbations and inflaton evolution 
are so significant, that care must be taken to account for these effects in the analysis of any inflation models. 

In recent work [J, |5| , we have identified a key mechanism which is generic in realistic inflation models and which 
leads to robust warm inflation. This mechanism involves the scalar inflaton field (f> exciting a heavy bosonic field x 
which then decays to light fermions ipd La > 

4>^x^^a. (l.i) 

In dynamical terms, this mechanism is expressed in its simplest form by an interaction Lagrangian density for the 
coupling of the inflaton field to the other fields of the form 
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£/ = -^9 2 <t> 2 X 2 - g'Hx^x " h X^d^d, (1.2) 

where ipd are the light fermions to which ^-particles can decay, with m x > 2m^ d . Aside from the last term in Eq. 
(|1.2|) , these are the typical interactions commonly used in studies of reheating after inflation [ill Il2f . However a 
realistic inflation model often can also have additional interactions outside the inflaton sector, with the inclusion of 
the light fermions ipd as depicted above being a viable option. Moreover in minimal supersymmetry (SUSY) extensions 
of the typical reheating model or multifield inflation models, the interactions of the form as given in Eq. i|1.2|) can 
emerge as an automatic consequence of the supersymmetric structure of the model. Since in the moderate to strong 
perturbative regime, reheating and multifield inflation models will require SUSY for controlling radiative corrections, 
Eq. 1|I.2[1 with inclusion of the i^d field thus is a toy model representative of many realistic situations. 

In Ref. [a] we have presented some results for the above mechanism in a fully expanding space-time dissipative 
quantum field theory formalism. The primary purpose of this paper is to supply the full details of the formalism 
used in Ref. |jg. This paper presents the various approximations used in the derivation, elaborates on the different 
aspects of dissipative evolution and radiation production and discusses the impact of these results on the density 
perturbation problem of inflation. Although the formalism presented here has general application, for most of this 
paper we will focus on the mechanism of Eq. I|1.I|) . The paper is organized as follows. In Sec. [H]we compare the 
inflationary dissipative process associated with Eq. (|1.1() to the reheating process after inflation in the standard cold 
inflation models. A simple picture of dissipation for both processes is presented. This should help to illustrate the 
differences between the two processes at work in the two cases. In Sec. IHII the Lagrangian model studied in this paper 
is presented. Also this section presents the real time matrix of Green's functions needed for our calculations and 
their evolution equations in expanding space-time. In Sec. II VI we present approximate solutions for these evolution 
equations using the WKB ansatz. The derivation of the effective equation of motion for the inflaton field from a 
response theory approach is done in Sec. [V] where we also discuss several simplifying assumptions for the nonlocal 
terms appearing in the effective equation of motion. These basic equations are then studied numerically in Sec. IVII 
and the validity of the different approximation schemes are explicitly tested for values of parameters of interest to 
inflation. In Sec. IVIII we discuss SUSY models that realize the basic interaction structure studied in this paper 
Eq. I|I.2|) . In Sec. IVIIII we study the effects this dissipative process has on density perturbations during inflation. 
Our concluding remarks are given in Sec. IIXI An Appendix is also included where some details are given on the 
renormalization of the effective equation of motion derived in this paper. 

II. INTERPRETATION OF DISSIPATION IN THE LINEAR AND NONLINEAR REGIMES 

It is useful to contrast the dissipation process to be discussed in the following sections to the dissipation process 
in the (old) reheating studies. For this, consider for instance the typical models for reheating, where an inflaton 
field <f>, with potential V (</>), is coupled either to spinor fields ip through the usual Yukawa coupling hfiipip and/or to 
other scalar fields x? with coupling g 2 (f> 2 x 2 ^ where in this last case the inflaton potential has symmetry breaking, with 
minimum at <p — 4>v- 

The typical reheating scenario is pictured in the time period just after the inflationary regime, where the inflaton 
energy density is released in the form of decay products of ip and/or \ particles. In the reheating regime the Hubble 
constant H is smaller than the inflaton mass m^, which means the inflaton can oscillate about the minimum, </>„, of 
its potential V(4>). In addition, to have particle production it requires m^ to be sufficiently large so 4> can decay, 
rrirf, > min(2m^, 2m x ). Typically one takes m^ 3> m^,, m x . In this regime the equation of motion is simple to treat for 
small inflaton amplitude. For example in a quartic inflaton potential with self-coupling A, a perturbative treatment 
is possible under the conditions A0(O) <C m^, m x , m^. In such a regime, the equation of motion for the homogeneous 
inflaton field </>(£), including quantum corrections, is given by the general linearized form 11] 

4>(t) + 3HJ>(t) + [m% + E(fc)] <t>{t) = , (2.1) 

where E(A:) is the polarization, or self-energy operator for </>, with four-momentum k = (oj, 0,0,0), with u> — m^. 
Due to the condition m^ > min(2m,/,, 2m x ) the <j) self-energy has a nonzero imaginary part ImE. The real part of 
the self-energy only renormalizes the mass m^, giving an effective mass to the inflaton, while the imaginary part is 
associated with the damping of <f) modes due to decay, with the decay rate given by T — — ImE(w)/(2u;). In the 
regime ImE <C m| , the solution of (|2.I|) ends up being the same as if we just replaced this equation with one having 
a friction like term proportional to T, 



(j>(t) + 3H<p(t) + ml4>{t) + T(f>(t) = . (2.2) 

Note that in this derivation, since we are considering the regime H <C m^, m Xl m^, the curvature of the universe is not 
important in the calculation of the self-energy. As such the decay rate calculation is typically just done in Minkowski 
space-time. Thus for the two decay processes of interest <f> — > x + Xi or — > "0 + V*) m the rest frame of the (^-particle 
the respective decay rates are 
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where we have considered m^ ^> m^, m x . Eq. I|2.2|l with Eq, I|2.3|) comprise the basic particle production process in 
the old reheating studies. 
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FIG. 1: The cut in the one-loop self-energy diagram and the amplitude of the decay process associated to it. Full lines stand 
for the scalar <j> and dotted lines to the spinors ip, tp. Analogous process follow for the decay <j> — > \X- 

As is well known [T3, the rates T also can be expressed in terms of the amplitude square for decay processes. This 
is illustrated in Fig. ^ where the imaginary term contributing to the rate of decays can be obtained by cutting the 
diagrams in half. In a dynamical problem, when the 4> modes of given energy 10 and momentum k are displaced from 
equilibrium, these modes get damped and it is the rate T that describes the approach to equilibrium. This approach 
to equilibrium is then naturally associated to a dissipative, or irreversible process. This picture is therefore closely 
related to that of a system (e.g. <fi) in interaction with an environment (e.g. the \ bosons and ip fermions), in which 
dissipation (and also noise) results from these interactions. More formally stated in quantum mechanical terms 14], 
we describe the state of the system by its reduced density matrix (by integrating out the bath degrees). Due to 
interactions, the state of the system gets entangled with the state of the bath and therefore some initial pure state of 
the system will end up turning into a mixture. This is an irreversible process, often called decoherence, that results 
from the non-unitary system evolution. 

We should also point out that although dissipation appears as a generic consequence for a system in interaction 
with a bath or environment, the simple representation for the dissipation, as in Eq. I|2.2|l and the derivation of it, 
implies a number of simplifying assumptions whose validity need to be checked. For instance, the derivation of Eq. 
(|2.2|) refers just to a very particular regime for the inflaton field, when it is oscillating around the minimum of its 
potential, with small field amplitudes and in the perturbative regime (or the linear relaxation regime). For this case 
the effective equation of motion of 4> has the simplified linearized form Eq. (|2.2|l . This simplified equation would 
not apply in nonlinear regimes, when large field amplitudes dominate the dynamics, for example in the description 
for the field modes during preheating, or in any other situation where nonperturbative effects play a relevant role in 
the description of field dynamics. Even in the linear regimes, Eq. I|2.2[| can be shown to be valid only up to a time 
interval At < 1/r [l5J |. beyond which the decay of cf> is no longer exponential but power law, which itself indicates 
the break down of the perturbative approximation used to derive Eq. (|2.2|) . It is clear that in more general cases of 
large field amplitudes or beyond the perturbative approximation, the expected effective equation of motion for the 
scalar field <f> must become very different than the simple equation (|2.2|l . Indeed, in general the effective equation of 
motion for an arbitrary scalar background field is a nonlocal equation. For example, in the nonlinear regime or for 
high field amplitudes, the description of the effective dynamics is not a simple local equation of motion (see e.g. Refs. 
). In addition, as we move away from the regime of validity of linear relaxation dynamics, it may become 
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possible to find other dissipative mechanisms not directly associated to the direct decay process that leads to Eq. 
(|2.2|) . In fact, as shown in recent work 0,0], even in the case where the inflaton can not decay, but fields coupled to 
it do (actually, <f> does not even need to be the heaviest field), dissipative regimes arise that are not available in the 
linear or perturbative regime. 

An example of interactions leading to a nontrivial inflaton dissipative dynamics is for instance the ones shown in 
Eq. (|1.2J) . where now m x > mm.(2m^ dl m ( p) and m$ < min(m^, ,ti x ). Here the \ particles can decay into fermions 
ipd that are coupled to it, but there are no kinematically allowed direct decays of <p into other particles. Nevertheless 
it is simple to understand from elementary particle physics the origin and nature of dissipation for the inflaton in 
this case. We look at those processes involving <p that may have an imaginary term and so in analogy to p.lfl . can 
be associated to dissipation. This is better interpreted in terms of an effective theory for <p after integrating over the 
other fields. We can start doing this by first integrating over the ipd fermion field. Since this field only couples to 
X, its main effect is to dress the \ scalar propagator, as shown schematically in Fig. [3 Note that the lowest order 
correction to \ goes exactly like the previous case analyzed above, in which the </> field could decay into ip particles. 
The leading order one- loop self-energy contribution to x, E x (fc), has a real part that represents a shift in the mass 
squared of the x and an imaginary part that represents the rate of its decay, as is kinematically allowed. This is the 
same kind of process shown in Fig. ^ by replacing the external lines by \ an d the internal ones by the spinors ipd- 
Next we can now perform the integration over the (dressed) x scalar particles and the ip x spinors coupled to (p. The 
relevant contributions to our dissipative mechanism are due to the x decay (see, however, later in Sec. V a discussion 
on the role of the ip x spinors). We now have processes contributing to the effective action to <p like the ones shown in 
Fig. |21 At leading order, the important contribution to dissipation in the effective equation of motion for <p> arises now 
from the one-loop vertex diagram shown in Fig. 0] By cutting that diagram in half we are now led to an imaginary 
contribution that can be seen as a dissipative term appearing in the effective equation of motion of (p. This can be 
understood from the amplitude shown in Fig. 0] which represents a scattering process of a by a virtual x that then 
decays into the fermion particles. This process can be interpreted in terms of the effective theory for (p, where an 
evolving background 0- field configuration excites </>-energy modes which then decay into the light fermions ipd, with 
that decay mediated by (virtual) x particles. The resulting dissipative term appearing in the effective equation of 
motion for (p can be seen from the square amplitude shown in Fig. 0]to be of order 0((p 2 g 4 T ^^^ . This result is 
in fact corroborated by the explicit derivation of the dissipation term in Sec. V. We can also note that this result 
is nonlinear in the <p field amplitude, since it originates from a scattering process (involving two <p particles and 
two virtual x particles) and it is nonperturbative in nature, since it involves the dressed x propagators based on a 
resummation, which means that physically the <p> field is not interacting with vacuum like x excitations but rather 
with the collective x excitations. 
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FIG. 2: The lowest order terms of the effective propagator for the scalar \ after integrating over the fermion fields coupled to 
it. The black ellipse stand for the spinors quantum corrections and the dashed lines to the x propagator. 

It should be noted that the dissipative processes as represented in Fig. 0|are of higher order than the ones shown in 
Fig. ^ However, due to the nontrivial nature of the latter, they may become important in those regimes characterized 
by high <p amplitudes (and then outside the region of validity of linear relaxation theory), and couplings (g and h in 
Eq. (|1.2|0 that are not perturbatively small. It is exactly in this region of parameters, that we will find the relevance 
to inflationary dynamics of the nonlinear dissipative mechanism discussed here and derived in this paper. 

III. A WORKING MODEL OF A SCALAR FIELD IN INTERACTION WITH MULTIPLE OTHER 

FIELDS 

Consider the following model initially presented in |i| , which consists of a scalar field $ interacting with a set of 
scalar fields %j, j = 1, . . . , N x and these scalar fields in turn interact with fermion fields ipk> k = 1, . . . , N$. Here we 
work in the FRW background metric ds 2 = di 2 — a(t) 2 dx 2 . The Lagrangian density for the matter fields coupled to 
the gravitational field tensor g^ is given by 





FIG. 3: Some of the lowest order diagrams contributing to the action of <j> after integrating over the scalar \ an d spinors ipd, ipd- 
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FIG. 4: The cut in the lowest order nontrivial diagram appearing in the effective action of (j> that has an imaginary part and 
the main amplitude associated to it. 
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where R is the curvature scalar, 
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(3.1) 



i? = 6 



(3.2) 



£ is the dimensionless parameter describing the coupling of the matter fields to the gravitational background and all 
coupling constants are positive: A, fj, g? hkj, x > 0. In the last term involving the fermion fields, the 7 M matrices are 
related to the vierbein e° (where g^ v = e®e^r] a b, with r] a b the usual Minkowski metric tensor) by [l]} 7 M (^) = l a e-a( x )i 
where j a are the usual Dirac matrices and w p = — {i/4)<j ah e'£V ll ebv, with a ab = i[7 a ,7 & ]/2. 

We are interested in obtaining the effective equation of motion (EOM) for a scalar field configuration (p = ($) 
after integrating out the $ fluctuations, the scalars Xj an d spinors ipk,tpk- This is a typical "system-environment" 
decomposition of the problem in which Lp is regarded as the system field and everything else is the environment, which 
in particular means the <£ fluctuation modes, the scalars Xj an d the spinors ip^, ijn. are regarded as the environment 
bath. In a Minkowski background, at T = 0, the EOM for ip has been derived in [3] using the Schwinger closed time 
path formalism. Here we follow a completely analogous approach and derive the EOM in a FRW background. The 
field equation for $ can be readily obtained from Eq. H3.lt and it is given by 



$ + 3-4 T $ + m\<$> + -$ 3 + £i?$ + 

a a z v 6 



5>f$(z)x^) 
.7=1 



0. 



(3.3) 



In order to obtain the effective EOM for ip, we use the tadpole method. In this method we split $ in Eq. I|3.3|l . as usual, 
into the (homogeneous) classical expectation value ip(t) = ($) and a quantum fluctuation <fr(x), <5>(x) = (p(t) + 4>(x). 
This way, the field equation for <f>, after taking the average (with (<f>(x)) = 0), becomes 



(p(t) 



3^|r0(t)+mfc(t) + ^ 3 (t) + 



ZR(t)ip{t) + ±y(t){ct> 2 ) + ±< 



Ny 



E^W<x?) + (^)]=° 



(3.4) 



3 = 1 



where (4> 2 ), (0 3 ), (x 2 ) and (4>Xj) can be expressed Q in terms of the coincidence limit of the (causal) two-point Green's 
functions Gt + (x,x') and G+ + (x,x'), for the $ and Xj fields respectively. These Green's functions are appropriately 
defined in the context of the Schwinger closed time path (CTP) formalism 18] . They are obtained from the (1, 1)- 
component of the real time matrix of effective propagators, which satisfy the appropriate Schwinger-Dyson equations 
(see, e.g., pHHI for further details). These equations satisfied by the effective (or dressed) propagators emerge from 
the successive integrations over the bath fields in l|3.1|l . The integration over the spinors ipki'ftk leads to dressed 
propagators for the Xj fields (see e.g. Fig. [3| , which are then given by (in the FRW background) 
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G x . (x, x') + / d 4 zS Xj (x, z)G Xj (z, x') = i 



8(x, x') 



a */*(t)a*/*(f) ' 



(3.5) 



where S Xi (x, y) is the self-energy for \j due to the coupling to the spinors ipk, ipk- Next, by integrating over the Xj 
and 4> fluctuations we are left with an effective propagator for the $ fields that is also formally defined by 



dt 2 a dt a 2 



$*t 



(R(t) 



G c j,(x, x) + / d 4 zT. ( i,(x 1 z)G<f,(z, x') = i 



5{x,x') 



a 3 / 2 (i)a 3 / 2 (i') 



(3.6) 



where Yl^(x,y) denotes now the self-energy for $ after integrating over the remaining bath fields. 

In the CTP formalism of quantum field theory Eqs. (|3.5(l and l|3.6|) are matrix equations for the propagators and 
self-energies, where, for instance, by expressing G^,(x,x / ) in terms of its momentum-space Fourier transform (and 
analogously for G Xj ), it can be expressed in the form 



G$(x,x') = i 



d 3 g 
(2tt) 3 



iq.(x-x' 



G++(q,t,t') G+-(q,i,t') 
G7+(q,i,i') G7-(q,t,f) 



(3.7) 



where 



G++(q, t, t 1 ) = G>(q, t, t')0(t - + G< (q, t, t')6{t' - t), 
G--(q, t, t') - G> (q, t, t')6(t' -t) + G< (q, t, t')6{t - t'), 
G+-(q,t,t') = G<(q,t,t'), 
G-+(q,M') = G>(q,M'), 



(3.8) 



with each matrix element defined in terms of two-point correlations of the fields that are on each branch of the CTP 
contour [l8J |. In writing Ij3.7|l we are assuming that G(x, x') only depends on the difference x — x', which follows for 
homogeneous field configurations, as is our interest here. The elements of the propagator matrix, Eq. I|3.8|l . for both 
the scalars bosons <fi and Xj, ap e found to satisfy the conditions 



G > (x,x') = G < (x',x) , 
iG >{< \x,x')Y =iG <(> \x,x') , 

G>«)(q,t ! t')] t = G>«)(q,t',t) ! 

-^[G>(q,<,0-G>(q,t',i)] =iS(t-t'), 

at l J t=t' 



(3.9) 



where the third condition in l|3.9[) is just a result of the first two conditions. The definitions of the retarded and 
advanced propagators are also given in terms of the matrix elements of the two-point function in the CTP formalism: 



G rct (x, x') = 6{t - t') [G>{x, x') - G< (x, x')] = G ++ (x, x') - G + ~ (x, x') , 
G adv (x,x') =6{t'-t) [G < (x,x')^G > (x,x')] = G ++ (x,x') - G~ + {x,x') . 



(3.10) 
(3.11) 



In particular, we have the known result that G rct (x,x') = G adv (x', x). 

The self-energy matrix elements are also expressed in a similar way to the propagators in the CTP formalism as 



££+(&, a/) - T>{x,x')6{t - t 1 ) + V<{x,x')e{t' - t), 
X~(x,x') = T>{x,x')6{t' -t) + V<{x,x')6{t - t'), 
S^ (x,x ) = — S^(a;, x ), 
S^ (x, x ) = — S^ (x, x ) , 



with analogous expressions for the \j self-energy elements. From (|3.12ll the following property follows, 



S++(.t, x') + S + -(x, x') + TT + (x, x') 



~(x,x') = 



(3.12) 



(3.13) 



In addition the elements in (|3.12|) are found in general to satisfy conditions similar to those in (|3.9|) . which are valid 
for both the scalar bosons <p and Xh 



T,>(x,x') =T, < {x',x) , 
jS>«)(a;,a;0l =iZ<>\x,x') 



(3.14) 



In terms of Eq. (|3.12l) and the equations for the fluctuation field modes derived from Eq. 1)3. lfl . we can write general 
expressions for the solutions for the G > ( < > propagator functions. For instance, consider the fluctuations equation for 
the Xj fields that can be obtained after integrating over the fermion ipk,tpk- This equation is obtained from the 
quadratic action in the Xj scalar fields, S^Xj] obtained from Eq. 1)3. 1|) . In the CTP formalism it is obtained by 
identifying fields in each branch of the Schwinger's closed time path contour, with fields in the forward and backward 
segments of the CTP time contour identified as xt ancl X7: respectively (see e.g. [l6J). In term of these fields, we 
can express the quadratic action for the Xj fields, after integration over the fermions as 



Si[Xj,Xj 



where M 2 



d 4 xa d (t) 




V 2 
a 2 (i) 



+ K 3 )xl 



dt 



d*xa A {t) J d*x'a 6 (t') - [ X f(x)^+;(x,x')xt(x')+xi(x)^-(x,x') X -(x') 

+ X ;(x)Z-+(x,x>)x+(x')+x;(x)V--(x 1 x')x-(x') 



V 2 
" a 2 (t) 



+ K x] 



(3.15) 



g 2 ip 2 (t) + £R(t), (p(i) is the background <& field, and T, Xj [x, x') denotes the fermion loop contri- 



butions, which dress the Xj fields. It is now useful to use in Eq. (|3.15(l the redefined fields given by [Tj| 



xf 



1 / + 

- 2 (xt+: 
z xj-xJ 



(3.16) 



along with the identity (|3.13|) , which lead to the result 



Si[x C j,X 



3' ^3 -I 



d 4 xa 3 (t) 



xt 



if 



d v 2 



TTTT -3~— + ^T -M 



<% 2 "oat 

d 4 xa 3 (t) I d 4 x'a 3 {t') {xf(x) 



Xj 



Kt 1 



■{x,x') + T,+7{x,x l ) - S-+(x,.x') - £—(*, aOj X,V) 
- xf (a:) [Sj+Cx,^) - S+-(x,.x') - £-+(*, *') + E—foaO] X f (a/)}(3.17) 
Using Eqs. (|3.12|) and (|3.13|) . the argument involving the self-energies in the second term in Eq. 1|3.17[) becomes 



^+{x,x') + T+-{x,x') 



Xj 



\x,x') 



Xj 



(x,x') 



= 2 



Xj 



(x,x') + T,+-{x,x') =26{t x -t 2 ) 



£>(x,a/) 



£<(*,*') 



n x . (x, x) = Ux tXj (x, x) + n 2 . X] {x, x) , 



(3.18) 



where 



IIi, X3 (x, x 1 ) = [26(h - t 2 ) - 1] [e> (x, x 1 ) - £< (x, x 1 ) 
n 2 , Xj (x, x 1 ) = £> (x, a;') - E< (a;, x') , 



(3.19) 



which, from (|3.14|l have the properties Hi tXj (x,x') = IIi ]X . (a;', x) and fl2. Xj (a;, a;') = — Il2. Xj (a;', a;). Similarly, the 
self-energy contributions in the third term in Eq. (|3.17(l can be written as 



£++(x,a;') -£+-(*, x') 



Xj 

2 



vjv-'-v ^x/C^^O + ^Xi ^'^ 



£^(a;,a;') + % Xi (x, x') 



2ilm 



Y,>.(x,x') + ^.(x,x') 



(3.20) 



where the last equality in 13.20|l follows from H3.14[) . from which we see that E?".(a;,a:') + E5-. (a;, a;') must be purely 
imaginary. 

By substituting Eqs. (|3.18|) and l|3.20|l in (|3.17|) . due to the result l|3.2U|l . we are led to an imaginary contribution 
to the effective action for the Xj fields. This imaginary term can be appropriately interpreted as coming from a 
functional integral over a stochastic field, which then turns the evolution equation for \j i n to a stochastic form 
due to the presence of a noise term jl(j . Thus the complete evolution equation for the modes and background fields 
includes noise terms. Though their study is particularly difficult, previous estimates of their effects on the background 
dynamics during inflation and reheating |2(| shows that changes in the dynamics and energy densities are marginal 
for chaotic inflation kind of models and within parameters values for coupling constants (X,gj, /ifc) corresponding to 
the cases of interest in this paper. Due to this we can neglect the stochastic noise terms appearing in (|3.17|) . On 
the other hand these noise terms are important in obtaining the first principles evolution equation for the fluctuating 
modes of the inflaton, such as for studying density fluctuations during inflation; however in this paper we will not go 
that far. Thus by defining the evolution equation for \j modes from 



g nit 
M 



5S 2 [X C 1 ,X?} 



= 



(3.21) 



we are led to the following equation for the \j modes f x - (q, t) in momentum space 



d 2 n ad q 2 J<r2/N 



f Xj (q, t) + J dt'a 3 (t')Il Xj (q; t, t')f Xj (q, t') 



(3.22) 



where Il Xj (q; t, t') is the spatial Fourier transform of the \j field self-energy term given by Eq. I|3.18|l . An analogous 

<i\ + ^ip(t) 2 +£R(t) and self-energies 



expression for the fluctuation <\> modes also follows, like Eq. (|3.22|) . with M?(t) = m 2 , 

terms coming from the dressing due to Xj an d </> loops. The initial conditions for these field mode differential equations 

will be explicitly stated below for the case of conformal time. Though this is not of special concern in this work, 



this is a convenient way to circumvent known subtle issues of renormalization dependence on the initial conditions in 
gravitational backgrounds when formulated in comoving time |21| . 

In terms of the general solutions of (|3.22(1 . /i,2(q, t) and their complex conjugate solutions, obtained equivalently 
from e.g. the complex conjugate of Eq. IJ3.22H (note that when the bath self-energy term entering in 13.22fl has an 
imaginary part, the equation becomes non-Hermitian) , we then can write general expressions for the CTP propagator 
terms, for both the scalars <j> and Xh m agreement with the continuity conditions expressed in Eq. Ij3.9|l . in the general 
form EllH 



G> (q, t,t') = h (q, t)h (q, t')6(t -t') + A* (q, t')/ 2 * (q, *W " *) > 
G < (q,t,t')=n(q,t)f^c l ,t')e(t-t') + h(ci,t')f2(c L ,t)e(t'-t). (3.23) 

The solutions /i,2(q, t) and the appropriated initial conditions needed to determine them are discussed below. 

IV. SOLVING FOR THE MODE FUNCTIONS AND REAL TIME INTERACTING PROPAGATORS 

Typically, equations for the mode functions for an interacting model, of the general form as given by Eq. H3.22JI . 
can be very difficult to solve analytically, in particular for an expanding background. There are, however, a few 
particular cases, like for de Sitter expansion H ~ constant, so a(t) = exp(Ht), and power law expansion a(t) ~ t n , 
where solutions for the mode equation for free fluctuations are known in exact analytical form (see e.g. Ref. |2l|). 
For instance the mode equation for free fluctuations in de Sitter is (where in this case the scalar of curvature becomes 
R = 12H 2 ) 



dt 2 



d_ 
di 



3H -^ + q 2 e- 2Ht + m 2 + 12£H 2 



/de Sitter (q, £) — , 

which has known solutions given in terms of Bessel functions of first and second kind, 



(4.1) 



/ 1;2 (q, t) ~ J v {qe~ Ht /H) , Y v (qe- Ht /H 



(4.2) 



with v = —i^Jm 2 /H 2 + 12£ — 9/4. The other case where we can find an exact solution for the modes corresponds to 
power law expansion, where, by considering a(t) = (t/to) n and massless (free) fields with minimal coupling (£ = 0), 
the solutions are given by 24] 



/i,2(q,*) 



t l/2 H il) 



qtfit 1 - 7 
n-1 



t l/2 H (2) 



qtfit 1 - 7 
n-1 



(4.3) 



with fj, = (1 — 3n)/[2(l — n)] and H^ ' (x) = J v (x) ± iY v {x) are the Hankel functions. 

Alternatively, for deriving an approximate solution for the mode functions in the interacting case, we can apply a 
WKB approximation for equations of the general form Eq. (|3.22|) and then check the validity of the approximation 
for the parameter and dynamical regime of interest to us. As will be seen below, under the dynamical conditions 
we are interested in studying in this paper, this approximation will suit our purposes. Let us briefly recall the 
WKB approximation and its general validity regime, when applied to obtaining approximate solutions for field mode 
equations. An approximated WKB solution for a mode equation like 



|r + - 2 (q,t) 



/(q,t) = 0, 



(4.4) 



, which holds under the general adiabatic condition 



is of the form fwKB(<l, t) = 1/ [w(q, £)] exp ±i J dt"uj(q,t") 

cj(q, t) <C w 2 (q, t). We must point out that there is no problems in extending this approximation to an expanding 
background and in fact it is a common approximation taken for instance in the analysis of perturbation modes in 
the adiabatic regime j25j. In that case, however, massless modes are considered and so the approximation holds only 
for large enough physical momenta q/a 3> H, corresponding to wavelengths deep inside the horizon. Here, instead, 



we work in the large mass scale regime, e.g. 



3> H. In this regime the WKB approximation is also valid. This 
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can easily be checked by comparing the WKB solution fwKB with the one obtained from the exact solution for the 
modes equations, e.g., given by l|4.2f> . f e xact(q,t), for free fields in de Sitter spacetime. For example, it is useful to 
examine the ratio I fwKFil fe ^rA- For both solutions the same initial/boundary conditions are taken (in conformity 
to the Bunch-Davies vacuum |21|) at to (and lets say to — 0), so both results match at the initial time. They also 
match in the asymptotic q — > oo or H — > limits, as they should, so as to correctly reproduce the Minkowski results. 
But they are also found to match very well for masses m ^$> H, independent of the value of the physical momenta (in 
particular even for q/a <-ff<m). For example, it can easily be checked that for (in units of H) £, — 0, q/a = Q.01H 
and for m = 10 H , the overall numerical difference between the exact and approximated WKB forms for the modes 
is at most not more than one percent for an evolution in the first 10 e-folds and this discrepancy gets smaller for 
longer evolutions. For the typical parameters we consider in this work (in Sec. VI) and relevant for the dissipation 

mechanism discussed in Sec. II, we have for example m x ~ 10 H, and so the WKB approximation is expected to be 
excellent, which is indeed confirmed by the numerical results to be shown later in this paper. In addition, note also 
that large mass scales, m ^> H, imply that curvature effects in the field quantum corrections to be considered for the 
background inflaton field are subleading, with the dominant terms being the Minkowski like corrections. 

Proceeding with our derivations, consider then a differential equation in the form of Eq. IJ3.22I1 . Instead of working 
in cosmic time, it is more convenient to work in conformal time r, defined by dr = dt/a(t), in which case the metric 
becomes conformally flat, 



ds 2 = a{r) 2 {dr 2 - dx 2 ) 
By also defining a rescaled mode field in conformal time by 



(4.5) 



a(r) 



/(q,r)=/(q,i), 



(4.6) 



we can then re-express Eq. (I3.22[) in the form (generically valid for either <fi or Xj scalar fluctuations) 



dr 2 



/(q, r) + uS(q, r) 2 /(q, r) + / dr'n(q, r, r')/(q, r') = , 



(4.7) 



where we have defined 



w(q,r) 5 



a{rf 



M 2 



£• 



G 



R(r) 



(4.8) 



In (|4.8|l the conformal symmetry appears in an explicit form, with £ = 1/6 referring to fields conformally coupled to 
the curvature, while £ = gives the minimally coupled case. Note also that in conformal time the scalar curvature 
becomes 



6 d 2 a 
R{T) = ^d^ ■ 

In Eq. (|4.7|l we have also defined the self-energy in conformal time as, 



(4.9) 



n(q, t,t') 

a{Ty/ 2 a{T>y/ 2 



n(q,t,0 



(4.10) 



where the self-energy contribution II, coming from the integration over the bath fields, is given by the space Fourier 
transformed form for Eq. (|3.18(l . In (|3.18() . II was split into symmetric and antisymmetric pieces with respect to its 
argument as defined in Eq. (|3.19|) Thus based on Eq. I|3.19fl . the self-energy term in (|4.7|) can then be written as 
ll(q, r, t') = IIi(q, t, t') + Il2(q, r, r'). In addition, by writing the self-energy term in a diagonal (local) form [261 1271 ] 



li(q, t,t') = n(q,r)<5(r - r') - [f^q, r) + li 2 (q, r)] 5(t - r') 



(4.11) 
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and from the properties satisfied by 111 and II2, it results that IIi(q,r) must be real, while D^q, t) must be purely 
imaginary. The real part of the self-energy contributes to both mass and wave function renormalization terms that 
can be taken into account by a proper redefinition of both the field and mass M. On the other hand, the imaginary 
term of the self-energy is associated with decaying processes, as discussed previously. So, we can now relate the decay 
width in terms of the CTP self-energy terms as 



Imll 

2d) 



2d) 



(4.12) 



and Eq. (|4.7|l can be put in the form 



-^ + d)(q, r) 2 - 2*w(q, r)f (q, r) 



/(q,r)=0 



(4.13) 



We can now proceed to obtain a standard WKB solution for Eq. (|4.13() . To do this, following the usual WKB 
procedure, we assume the solution to have the form /(q, r) = cexp[i7(q, r)], where c is some constant that can be 
fixed by the initial conditions, given by 14.18JI below. This form of the solution is then substituted into (j4.13|) to give 



h" - i 2 + u 2 - 2iu>r = 



(4.14) 



Working in the standard WKB approximation, for the zeroth order approximation we neglect the second derivative 
term in l|4.14Jl . Then, by taking T <C Q, we obtain 



7o « T / dr' (d) - if) , (4.15) 

J To 

which is then used back in (|4.14() for the second derivative term to determine the next order approximation, 



; 1 ~- = / dr' [uJ-iT + O (d)' 2 /d) 3 )] + i In \/d) . 
't 



(4.16) 



The next and following orders in the approximation brings higher powers and derivatives of w'/w 2 , which in the 
adiabatic regime, d)'/d) 2 <C 1, are negligible and we are then led to the result 



/i,2(q,r) « -^exp 



=pi / dr' (oj — ir) 

TO 



(4.17) 



The solutions for the modes of the form l|4.17|l and their complex conjugate are general within the adiabatic, or WKB, 
approximation regime of dynamics. Finally, we completely and uniquely determine the modes by fixing the initial 
conditions at some initial reference time r , which can be chosen such that in the limit of k — > 00 or H — > we 
reproduce the Minkowski results. These conditions, which correspond to the ones for the Bunch-Davis vacuum |21j . 
can be written as 



vM T o) 

/ i,2(q, T o) = tV^( t o)/2 , 



(4.18) 



which already fixes the constant c in (|4.17|l asc= l/v2- 

Using the above results in l|3.23|) and after returning to cosmic time t, we obtain the result, valid within the WKB 
approximation, or adiabatic regime, 



G>M(q,M') 



[a(t)a(t')} 



3/2 



G>«>(q,M') 



(4.19) 
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where 



G>( q ,M') - f ,/, ^ (e"*^ d '"[-('")-^(*"% (i _ i') + e ~*/; ««"[«(*")+ir(t")] < , (t / _ t) 

VM ; 2[w(£)w(t')] 7 I 

G<(q,t,t') = G>(q, £',£), (4.20) 

where T is the field decay width in cosmic time, obtained from (|4.12|) and 



«(*) = J ^ + M*(t) , (4.21) 



with M 2 (t), for $ particles, given by 



M}{t) = m\ + \ V (tf + (t~l) R(t) , (4-22) 

while for \j particles, 

Ml (t) = ml + 3 >(t) 2 + ($-") ^(*) • (4-23) 

The same result (I4.20f) could in principle be inferred in an alternative way by expressing the propagator expressions 
in terms of a spectral function, defined by a Fourier transform for the difference between the retarded and advanced 
dressed propagators, given by Eqs. (|3.1UI) and i|3.11f) . and approximating the spectral function as a standard Brcit- 
Wigner form with width given by T and poles determining the arguments of the exponential in (|4.17() and its complex 
conjugate [lq . The validity of this approximation in particular was recently numerically tested and verified in Ref. 
[28| for a 1 + 1 d scalar field in Minkowski space-time. In the Minkowski space-time case, results analogous to Eq. 
(|4.2U|I were explicitly derived in Refs. |ll l3L ll6L Il9| . Indeed, for the case of no expansion a(t) = constant, Eq. (|4.2U|) 
reproduce the same analogous expressions as found in the case of Minkowski space-time. 

The result <|4.2U[1 . from the previous approximations used to derive the WKB solution 14.17fl . is valid under the 
requirements 

r^ < ^0 , 

T Xj « Uxi , (4.24) 



and the adiabatic conditions, 



■jj 



_a/a | L)^ 1 



w 2 u<t> ujI 



w v - a/a Cj-,, , 

-^ = — J- -g- < 1 , 4.25 



^ 7 ~\ . W 2 



Xi x ' X 3 

where in the second term in the equations Q4.25J I we have made the change back to comoving time and used uo — 
a(t)u)(t) = a\J q 2 ja 2 + M 2 . The conditions 14.2411 are generically valid in perturbation theory. The second set of 
conditions given by l|4.25[) are the usual conditions imposed in the derivative expansion for the WKB solution (which 
recall was here obtained for convenience in conformal time). These conditions are valid whenever the adiabatic 
conditions for the background field ip(i) are satisfied, which is the case for a slowly moving field. They are also 
satisfied for those modes deep inside the horizon, q 3> aH, which is useful when expressing the WKB solution as a 
large momentum expansion and for explicit renormalization purposes. Finally, the condition Q4.25J I is also found to 
be satisfied for those modes outside the horizon, q <C aH, provided the masses M$ and M x . are much larger than 
the Hubble scale and their time dependence evolves in an adiabatic manner. To better sec these different regimes of 
validity of Eq. 1)4. 25|) . we look at the two extreme cases of parameter regimes of interest in this paper. For both these 
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cases M ^> H . In the first extreme, the modes are outside the horizon, q/a <C H <C M, in which case Eq. (|4.25(l 
becomes 



H£ ~ T^ « 1 



d/a cJ0 

W0 W^ 


M 


w x, <*& 





— + -T- ~ TTT « 1 . ( 4 - 26 ) 

Uv,. W?. M. ^ ' 

Xj Xj Xj 

For the other extreme case, the modes are deep inside the horizon, q/a ^> M 3> H, in which case the conditions Eq. 
(I4.25|> automatically become satisfied, since 



UJ4, w\ 

a/a u Xj q/*>M X j>H 



. (4.27) 

u x 3 <• 

This last case is the weakest condition, since during inflation the tremendous growth of the scale factor makes the 
modes rapidly go outside the horizon, thus going over to the regime of the first set of constraints Eq. (|4.26|) . Condition 
Eq. H4.26f) can be satisfied provided the background field moves sufficiently slowly, which is the regime we will be 
interested in probing in this work. It should also be noted that in the parameter region of masses M 3> H, curvature 
effects become sub-leading and so Minkowski like expressions can apply to leading order. We therefore expect that in 
the adiabatic regime the approximations used to derive Eq. (|4.20() readily hold. This will be tested numerically later 
on in Sec. IV1I 

V. DERIVING THE EFFECTIVE EQUATION OF MOTION FOR THE INFLATON 

We now turn our attention to the EOM Eq. 1J3.4I) , where we will work it out in the response theory approximation 
similar to the treatment in our recent paper [4J . Consider the Lagrangian density (|3.1|1 in terms of the background 
(system) field ip(t) and the fluctuation (bath) fields, 

£[$ = ip{t) + c/)(x),Xj,-ipk,ipk,g^} =£<p[<p(t),g liV ] + C hi ,th[ip(t),(f>(x),Xj,^k,tpk,g^} , (5.1) 

where 



C v W), 9liV ] = a(tf -m 2 - -^fitf - - V {tf - \Rip{tf , (5.2) 





is the sector of the Lagrangian independent of the fluctuation bath fields, while £bath denotes the sector of the 
Lagrangian that depends on the bath fields. In the following derivation it will be assumed that the background field 
ip(t) is slowly varying, something that must be checked for self-consistency Thus, if we consider the decomposition 
of (p(t) around some arbitrary time to as ip(t) — <p(to) + ^f(t)> fip(t) can be regarded as a perturbation, for which a 
response theory approximation can be used for the derivation of the field averages in l|3.4|) . I n order to implement 
the response theory approximation, we consider the terms in £bath that contribute to the derivation of those field 
averages in the tp-F,OM Eq. (|3.4|) and take ip(t) = ip(t ) + 5<p(t). We denote those terms that depend on 5ip(t) as £ ir ^, 



~[2pfo)M*) + M*) 2 K-^' 



4Tt = «t)*{-$[Mto)6rtt) + 6<p(t) 2 ]<f?--£8<p(t)<l> 3 



1 5<p 



4 [MtoWv® + MO 2 ] x* - gpviWx 2 



(5.3) 



and we treat these terms in C^ t as additional (perturbative) interactions. 
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A. The Response Theory Approximation 



In response theory we express the change in the expectation value of some operator 0(t), S((D(t)) = (0(t)} peT t — 
(0(t)), under the influence of some external perturbation described by -ff Dcr t which is turned on at some time to, as 
(for an introductory account of response theory, see for instance Ref. |29j ) 



S(0(t)) =i / dt'( 

Jt 



H pert (t'),6(t) 



(5.4) 



where the expectation value on the RHS of i|5.4|) is evaluated in the unperturbed ensemble. The response function 
defined by Eq. I|5.4|) can be readily generalized for the derivation of the field averages. Provided that the amplitude 
5tp(t) is small relative to the background field (p(t ), perturbation theory through the response function can be used to 
deduce the expectation values of the fields that enter in the EOM Eq. i|3.4|) ■ In this case the perturbing Hamiltonian 
flpert is obtained from C^ t , Eq. Q5.3JI . From Eqs. Q5.3JI and l|5.4|l we can then determine the averages of the bath 
fields, for example ((j> 2 (t)), as an expansion in 5ip(t), starting from the time tg and in an one- loop approximation, as 



b 2 ) Q -i f dt'a{t'fj [2^(i )Mi')+<W) 2 ] ([0 2 (x,i),0 2 (x,£')])o + 0(<V 3 ) 

JtQ 



(5.5) 



where (. . .) means the correlation function evaluated for the background field taken at the initial time, (. . .) = 
(...) . The first term in (|5.5|) is just the leading order one-loop tadpole term in the linear response approximation, 

<p(to) 

while the second one is the one-loop tadpole made up with the interaction term from l|5.3|) . — a(t) 3 (A/4) [2ip(to)5ip(t) + 
5ifi(t) 2 }4> 2 , that is used in calculating the leading order one-loop bubble diagram that gives the two-point function. 
The interaction vertex coming from the above term can also be put in the more convenient form — ia(i) 3 (A/4) [ip(t) 2 — 
<p(to) 2 ], which we will use in evaluating l|5.5|l . 

Using translational invariance we can now write ([</> 2 (x, t), 2 (x, £')]), in Eq. I|5.5ll . in terms of the causal two-point 
Green's function for the <\> field, G\ + (x,x'), as 



([0 2 (x,£),0 2 (x,O]> = 2iIm(T^ 2 (x,^ 2 (x,0) 



= Ai 



rf 3 c 



a{t)a(t')Y J (27T) 3 l ^ 



Im[G+ + (q,i,t')] 2 >t , , 



(5.6) 



with Gt + (q, t — t') as obtained from Eqs. I|3.8() and l|4.20|l . Eq. I|5.5() in the response approximation then becomes 



a{tf 



dfx Up{t'y - tpfor 



to 



d 3 q 

(2tt) 3 



Im 



Gt + (ci, t,t') 



v(to) 



t>t' 



Analogously for the other field averages we find: 



and 



<X?) - (X 2 >0 + 



y " " V; " ' a(i) 3 



dt'2g 2 [ipit'f - V (t ) 2 



(</» X 2 ) ~ (0x-)o + 



a (ty- 



dt'2g 2 Mt') 



d 3 q 

J2^f 

1 



Illl 



G++(q,t,0 



J ^ K ' [a(t)a(t')] 3/2 J (27r) 3 (27r) 3 



v(to) 

d 3 qi d 3 q 2 



t>t' 



x Im 



G+ + (c Ll ,t,t')G++(< l2 ,t,t')G++(q 1 + q 2 ,t,t r ) 



t>t' 



P(*a) 



(5.7) 



(5.8) 



(5.9) 



<tf 



dt'2\<p(t'y 



d 3 qi rf 3 q 2 
a(i)a(i')] 3/2 J (2tt) 3 (2tt) 3 



x Im 



G+ + ( qi ,t,OG+ + (q2,i,t')G+ + ( qi +q 2 ,i,<') 



t>t' 



v(*o) 



(5.10) 
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Eqs. H5.7J I. (|5.8|l . (|5.9|l and (|5.10() are analogous to the ones obtained in 3] but derived there in the context of 
the Schwinger's closed-time-path formalism in Minkowski space. The leading order terms in the linear response 
approximation, (4> 2 }o, (~X 2 )o, etc, are divergent and need appropriate renormalization in expanding space-time, as 
e.g. described in Ref. |3tj; below we will give the explicit expressions for the relevant terms. While the first two 
expressions, Eqs. I|5.7() and l|5.8[) correspond, when expressed diagrammatically, to the one-loop tadpoles of one and 
two vertices in the ip-FiOM, the last two expressions, Eqs. (|5.9(l and (|5.10J) . represent two-loop contributions to the 
EOM. In the following, as in Ref. [j], we will restrict our study of the EOM at one-loop order. This makes our 
analysis tractable and simple. Moreover there is no loss in our analysis of the dissipative dynamics for cp, since the 
one-loop terms will already suffice to demonstrate the possible different dissipative regimes and the higher order terms 
only enhance the dissipation effects obtained in the analysis that follows. 

B. The (^-effective EOM 

We then obtain that the EOM Eq. (|3.4() . with bath field averages evaluated from the response function and at 
one-loop order, becomes 



^< t »^/(0 <5J+(q ' t ' < i<,» + l^ <() ^p/l?^ (q,M) 



v(*o) 



A . . 1 
-o ¥>(*)- 



1 T w a{tf 



ta 



dt' A [ip(t') 2 - v {t Q ) 2 



d 3 q 
(2tt) 3 



Im 



Gj + (q,t,i') 



-I 2 



■E^)^l o d^W) 2 -<P(t ? 



d 3 q 

J2nf 



Im 



v(*o) 
G++(q,t,f 



t>v 



¥>(*o) 



(5.11) 



t>t r 



We can now use Eqs. <|3.8[l and i|4.2L)[l and the equivalent expressions for the \j propagator in the above equation to 
obtain 



6 



d 3 q 1 



A 1 

'2^'^itf J (27r) 3 2^(i) 



<p(t ) 



sr ^ m i f d3 i l 



(t) 



¥>(*o) 



X 2 



N, 



<p(t) / dt' [rtt'f - <p(t ) 2 ] D^t, t') - ]T 2 »>(*) / dt' [<p{t') 2 - p(t ) 2 ] D Xj (t, t') = 0, (5.12) 



to J=1 



where the kernels D^(t,t') and D Xj (t,t') in the above equation are given, respectively, by 



D*(t,t') 



a(t) 3 J (2tt) 



d 3 q . 

■ sin 



2 / dt"io lP (t") 
t' 



exp 



-2^t%(q,f) 



4w^(t)w^(f) 



t>t' 



(5.13) 



and 



D x At,t') 



d 3 q 
a{ty J (2?r) 3 



2 / dt"LO x {t") 

v 



exp 



-2/ t !di"r x (q,t") 



4w x (t)w x (f) 



t>t' 



(5.14) 



where Y ^ and r Xj . are the decay rates for the <& and \j particles with momentum q, respectively. These depend 
explicitly on the decay channels available for both the $ and Xj fields, within the kinematically allowed masses. 
These we will fix explicitly below. Note also that, as a consequence of the linear response approximation, all the 
frequencies appearing in the above expressions are expressed in terms of ip = <£>(£o)> 



u<p(t) 



q 2 /a(t) 2 +ml + -^(i ) 2 + (£ - l/6)i?(i) 



1/2 



w« (t) = q 2 /a(0 2 + "& + fl V(*o) 2 + (£ - l/6)fl(i) 



1/2 
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(5.15) 
(5.16) 



In Sec. IV CI we show how the next order corrections in the linear response approximation (at one- loop order) can be 
resummed to give back the full time dependence for <p inside the above expressions. 

Eq. I|5.12J1 . with Eqs. Q5.13J I and Q5.14J I. is our general expression for the one- loop effective EOM for the background 
(inflaton) field ip. As in the Minkowski space case yl lJ we expect that the last two, nonlocal terms in Eq. (|5.12|) 
will lead to dissipation. This can be made apparent once we integrate them by parts with respect to t' . This way we 
obtain explicitly first order (nonlocal) time derivative terms in the background and separate additional local terms 
that, when combined with the first two momentum integral terms appearing in Eq. (|5.12ll . will correspond to the first 
derivative, d/dip, of the one- loop quantum correction to the effective potential V c ff(tp) (in the equation 1)5. 12|) we have 
the corrections from both the scalar field $ self-coupling and due to its coupling to the Xj fields). In the absence of 
the nonlocal dissipative terms and the additional couplings to <f>, this way of obtaining the (field derivative of the) 
one-loop effective potential was shown explicitly by Semenoff and Weiss in 22] and its renormalization later studied 
by Ringwald in |30|. 

To make more transparent the interpretation of the different terms that can be derived from Eq. I|5.12|l . let us 
define here dissipative kernels K^,(t,t') and K Xj (t,t') as related to the kernels Eqs. I|5.13[) and (|5.14|) . respectively, by 



dt> 



= D^(t,t'), 



(5.17) 



and 



dK Xj (t,f) 



dt' 



D x At,t'), 



(5.18) 



whose solutions we choose here so that in the limit of flat space and as we take in 



the kernels Ka, and K 



become the ones obtained in Minkowski space calculations 3] . From this we then have the solutions 



Xj 



and 



K (/> (t,t') = [ drD^r) , 

J t 



K Xj (t,t')= / drD Xj (t,r) 

Jt 



Using Eqs. I|5.17|l and I5.18fl in I5.12fl we obtain that 



(5.19) 



(5.20) 



<p(t) + 3^0(t) + m%tp<f) + ^(t) 3 + £R(tMt) 
a(t) v 6 



A J_ f d 3 q 1 

2 n, a(t) 3 J (2tt) 3 2 W0 (i) 



Wo) ^ 






1 f d 3 q 1 



Av 



a{tf J (27r) 3 2w x (t) 



v{to) 



JV V 



- ]T 2gfy(t) [pit? - p(t ) 2 ] K Xj (t, t) 



i=i 



-AV(*) / dt'v{t')p{t')K^t,t>)+Y. A 9 4 M t ) / dt'<p(t')<p(t')K Xj (t,t')=0 



(5.21) 
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C. The local terms and the effective potential corrections 



Note that the local terms in the second and third lines in Eq. I|5.21[) can be written as a field derivative of the one- 
loop quantum corrections from the Xj an d </> scalar field fluctuations, to the effective potential for the ip background 
configuration. This is easily seen by writing the local terms as 



<1 



1 



A i r d 3 

2 m ^tf J (2^)3 2w^(t) 



-— ip(t)[ip(t) 2 -ip(t ) 2 ]K^(t,t) 
tp(t ) i 



E 2 , . 1 f d q 1 
d 3 q 1 



JV V 



v(*o) 



X>^(t)[p(t) 2 -p(to) a ] #*(«,*) 



i=i 



2^*M a(i) 3 7 (2tt) 3 2^ 



-A [^(i) 2 - ip(t ) 2 } -i-3 
v(*o) a(t) 3 



(2tt) 



&4 



¥>(*o) 



+ °(i*/<4) 



+§**> {^/ 



d 3 <? 1 



a(i) 3 7 (2tt) 3 2cj Xj 



*,-^W ! » 2 -^]^?/|p 



W, 



f{ta) 



+ °<&J«>1) 



>5,22) 



where the C(r^/w^) and 0(r 2 ./u') explicit contributions would correspond to higher than one- loop contributions 



resulting from the use of the full (resummed) propagators obtained from Eqs. I|3.6(l and (|3.5|l . Eq. Q5.22J I can now 
easily be recognized as originating from the dip amplitude expansion of the local (free) propagators G^(t,t) and 
G x -(t,t), respectively, 



Gi + (t,t) 



( fc 



a{tf J (2tt) 



(fc 



q 2 /«W 2 + ml + ±p(ty + (£ - l/6)R(t) 



a(t) 3 J (2tt) 3 2w 



A [ip(t) 2 - p{t f 



»(*)* 



1/2 



d 3 q 1 

(2^R 



C((S^ 3 



(5.23) 



and 



GiUt,t) = 



i 



«W 3 
i 



d 3 c 



1 



(27T) 



A 



q 2 /a(i) 2 + ml , + ff 2 ^) 2 + « - l/6)fl(i) 



1/2 



1 



a(t) 3 7 (2tt) 3 2w : 



Xj 



- 2. 9 2 [^(i) 2 - ^(i ) 2 



1 



dh 



i 



a(i) 3 7 (27r) 3 8w| 3 



o(V) ■ 



(5.24) 



This confirms our above statement that these terms arise from the field derivative of the (unrenormalizcd) one-loop 
effective potential quantum corrections coming from the $ self-interaction and Xj coupling. Combining the result 
(|5.22|) with the tree level part of the potential and using (|5.23|l and l|5.24|l , we can then write 



dV cS (ip,R) 
dp 



JV V 



m%pit) + '-ip{tf + ^R(t)ip(t) + -ip(t)G++(t,t) + Y.9 2 Mt)G++{t,t) . 



6 



(5.25) 



i=i 



The two last terms in Eq. (|5.25(l are of course ultraviolet (UV) divergent as expected and so need to be properly 
renormalized. This is done in the usual way by adding to the original Lagrangian, or in the effective EOM for ip, the 
appropriate countertcrms of renormalization, dm^, SX and 5£, for the mass, scalar $ self-coupling and the gravitational 
coupling, respectively. The details of this renormalization process are discussed in the Appendix A, where explicit 
evaluation and renormalization are done. After renormalization we can just rename the couplings and masses in the 
dissipative and quantum corrections as the renormalized ones. Note also, as evident from Eq. I|5.25|l and the explicit 
results shown in Appendix A, that the quantum corrections can be kept relatively under control for perturbative small 
couplings (and small number of fields). This is certainly true for the contributions coming from the $ scalar field 
self-coupling, associated to the inflaton, which is required to be tiny (A < 10~ 13 ) due to the density perturbations 
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constraints requiring a very flat potential. However, for the x field contributions there are no physical constraints that 
require the couplings gj or the number of fields N x to be sufficiently small. In fact, in Refs. [3,|IJ relevant scenarios 
of strong dissipation are found for cases of intermediate to large couplings 0(1O~ 4 ) < N x gj < 0(1). In these cases 
we must worry about the large quantum corrections appearing in (|5.25|) and in particular that they will not spoil the 
flatness of the potential. We here follow the same procedure adopted in Refs. 0,01 to overcome this problem and add 
to the original Lagrangian density ((3.1(1 an additional coupling of $ to N$ extra fermion fields, — X}j=i di^^Phx^hX' 
where tp x are fermion fields, which are different from the light ones coupled to x i n ffi-lfl - For appropriately tuned 
coupling g' ~ g this modification just mimics supersymmetry, where $ couples to both the boson scalars x an d their 
fermion partners, with large cancellations occurring between the quantum corrections from the x an d V'x fields. This 
can be seen explicitly in V c r (</?), when the ip x fermion coupling to $ is included, which to one- loop order gives 

m l .2 , V.4 , £ 



V cS (f,R) = -^<p 2 + -y + \Ry 2 




N 



*i 






2 J (2tt) 3 '"* ^ '"* J (2tt) 3 ^ 

where fe p = ft/a is the physical momentum and 



£ mj — 



^/k p 2 +m2+A^/2+(£- 1/6)jR> 



£„, = Jkp 2 + m a + . 9 2 ^ 2 + K - 1/6)^, 



E 1 .,, 



, 'i/j i 



v /k p 2 + (m^, x+3 »2. (5.27) 



Thus, with appropriately tuned parameters gi, ^ and with zero explicit masses m t f, i — m Xi — and N^ i = N x /4 1 , 
the one-loop quantum corrections to V e s cancel to all orders in g t , gf in the nonexpanding case (R = 0). Even when 
a vacuum energy is considered (e.g. in de Sitter where R = 12iJ 2 ), the combined nonvanishing contributions from % 
and ip x can still be made small enough compared to the three level potential [3l]]. Further discussion about SUSY 
models is in Sec. VII. 

Note also that adding the interaction term — J^i=i 9fe' l l , i,x' l l J i,x (with the corresponding kinetic term for the ad- 
ditional fermion species) to the original Lagrangian density 1)3. l|l will also yield additional contribution to the EOM 

for ip, given by Y)j—x 9i{' t Pi,x' t Pi,x)- This term can be worked out analogous to the scalar case. Aside from the local 
corrections discussed above, this term will lead to an additional dissipative kernel in l|5.21|l . As shown explicitly in 
J3J, where this term was derived, it will not be directly proportional to the <p field amplitude, unlike the two nonlocal 
terms in (15.21(1 coming from the scalar field quantum corrections that are directly proportional to the square of ip. As 
in the old reheating like scenario discussed in Sect. [HI the fermionic nonlocal term will then be relevant in the linear 
regime (or small ip amplitude), while the last two terms in (|5.21() will contribute mainly in the nonlinear regime of 
interest here. Thus, we can just restrict our following analysis of the dissipative kernels to the ones given in ((5.21(1 
and neglect the contribution to the dynamics coming from the ipi% interaction, keeping in mind that the addition of 
any other bath fields coupled to $ will also add to dissipation or nontrivial effects that can play a role in different 
dynamical regimes. We will briefly return to this again in the conclusions in connection to the description of late time 
effects in the dynamics of p. 

D. The effective nonlocal EOM and energy densities system of equations 

With the considerations above, we can now write the effective equation of motion for the background field ip as 



1 When the \ and ip x couplings to <I> are treated in a SUSY context, this last restriction on the number of fields is not important, since 
X are then actually complex fields while their fermionic partners are Majorana spinors and the two contributions appear in 15.261 with 
the same number of degrees of freedom. 
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<p(t) + 3H(t)ip(t) 



dV c r s (<p(t),R(t)) 
d(p(t) 



iVv 



ft "x ft 

+AV(t) / dt' <p(t')<p(t')Kt(t,t') +Y t *8*M*) / dt'<f(t')<f(t')K Xj (t,t') = , 

J to j — l J t 



(5.28) 



where V^g stand for the renormalized effective potential (see Appendix A) and, from l|5.13[l . (|5.14|) . (|5.19|l and 15.20f> . 



K4t,t')= / dr 



it, 



d 3 q 



a{tf J (2tt) 



■ sin 



dt"u><t>(t") 



exp 



-2/>"I>(q,i") 



and 



K X] (t,t') 



dr 



a{tf 



d A q . 
■ sin 



(2tt) 3 



di w x [t ' 



exp 



4w0(i)w^,(r) 



-2^dfT x (q,t") 



v>v 



(5.29) 



4w x (t)w x (f) 



t>t' 



(5.30) 



are the dissipative nonlocal kernels. 

The complete evolution of the inflaton field is then determined from Eq. (|5.28|l and the Einstein equations for the 
background cosmology. Together, these equations form a complete set of dynamical equations for both <p and the 
metric. The Einstein equations for the background cosmology can be formed in terms of the matter and radiation 
components as usual by the equations: 



H 2 = 



iwG 



(Pm + Pr) J 



(5.31) 



and 



2H + 3H 2 



-8irG(p m + p r 



(5.32) 



where G — 1/mpj, with mpi the Planck mass. The parameter k = 0, +1,— 1 for a flat, closed or open Universe, 
respectively. In this work we only consider the flat case, k = 0. p m M and p m / r ) are the energy and pressure densities 
for matter (radiation), respectively. We also have the standard relations: 



tf + V&faR), 



(5.33) 



Pm = -if 2 -V^{^p,R) 

and p r — |p r . 

The matter and radiation energy densities p m and p r evolve in time as: 



(5.34) 



Pm + 3Htp 2 +\ 2 <p(t)<p(t) / dt' (p(t')(p(t')K^t,t') 



"x ft 

^4gj<p(t)<p(t) / dt^(t')f(t r )K X] (t,t') = 

3=1 *° 



(5.35) 



and (from the energy conservation law) 



v2. 
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p r + 4:Hp r ~\ 2 tp(t)<p(t) d£ ip{t')ip(t')K<i,{t,t') 

Jta 

- J2 4ff>(*)vK*) / dt' v (t')v{t')K X] (t, t')=0. (5.36) 

Assuming a flat universe (fc = 0), from Eq. I|5.31[l . we can consider 

pr = ^g r2 ~ pm = ^g h2 ~^~ v:s{ ^ R) (5 ' 37) 

as the first integral of Eq. (|5.36|l . Using Eqs. (|5.31() - (|5.34() . we can also express the equation for the acceleration in 
the following form 

- = ^-( Pm ~Pr)-^Gv 2 . (5.38) 

a 6 

The above equations together with Eq. (|5.21[1 form a closed, general set of (integro-) differential equations for 
the effective evolution for the background field ip(t) and metric at one-loop and leading order in the linear response 
approach. 

E. The equation of motion in a local approximation 

The derived equation of motion for (p, Eq. (|5.28|1 . in expanding space-time is a considerably more complicated 
expression than e.g. the analogous one that would be derived in the Minkowski case. It is therefore interesting first 
to see when and whether we can recover expressions equivalent to the Minkowski space ones in Q. This can be 
the case for instance if we restrict the dynamics in the adiabatic regime close to equilibrium, which then requires in 
general that the decay rates are larger than the Hubble constant, T ^S> H. Using this and restricting to time intervals 
t — t' ~ 1/I\ where the scale factor consequently changes very little, the frequency terms in the inner time integrals 
inside the kernels expressions Eqs. I|5.13|l and (|5.14|) will change very little, and they can be taken just as constant 
terms. Under these circumstances we can then easily see that Eqs. I|5.13|l and (|5.14|) can be approximated to 

D ^ t] ^tf] J^^l^m-t\] ^y 2 , (5-39) 

and 

r. / n 1 f d 3 q • r„ /m in exa[-2T y (Q,t)\t-f\] , „„. 

d <m~wJ ( 2 ^ sin ^ (t)|t - i| ] 4V — ' ( } 

which are equivalent to the kernels derived in 0. In terms of l|5.19|) . [f5~2T)jl . (|5~39"|) and JOGJ, the EOM Eq. (|5~2T|) 
now becomes 

m + aff^t) + dv ^ R) 

dtp 
2 ,n(A C rlt< <n(t'\,h(+>\— f ^ q t* cos(2^\t - t'\) + I> sin(2^|t - t'\)} 2r Jt _ t ,, 



+X 2 ip(t) / dt' ip(t')tp(t')- 



to 



a{tfj (2,)3 ^(rj+uj) 



e 



N- 



X- 4 M /"*.,/ ( ,s. f ,, 1 f d 3 q [wx, cos(2^ Xj \t - t'\) + r x . sin(2a; x . \t - t'\)] _ 2 r x .| t _t'] . ,. , n 
y,9jP{t) dt ip(t )<p(t — -J / TTTTE 7 — \ e J '=0 5.41) 

U J ^ «W 3 i(2-) 3 2<.(r X3+< ) 
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In this paper we restrict the analysis to zero temperature (or in a non-thermalized bath) and use the same regime of 
parameters as studied recently in |4j, where the masses (the renormalized and, if relevant, background field dependent 
ones) satisfy the condition M Xj > 2M^ k > M^. In this case T^ vanishes, while T Xj gives the decay width for the 
kincmatically available decay channel of the scalar \j fields into the fermion fields ipk,tpk- To obtain an expression 
for T Xj , first note that under the condition M Xj 3> H, which is satisfied for the couplings and values of the inflaton 
amplitude taken here, the curvature effects become negligible in the computation of T Xj and, therefore, it can be well 
approximated by the Minkowski decay rate at leading order. Thus, we can write its expression in terms of the decay 
rate in the rest frame, T x . (0), and then boost it to give (a similar form for the decay of a scalar particle into fermions 
in de Sitter space-time was also used by Ringwald in Ref. 26]) 

r xi(*) = ^% r x,(0)> (5-42) 

where T Xj (0) is the standard on-shell rate, as evaluated in Minkowski space-time J3J , 

N * M ( AM 2 \ 3/2 

This result quoted for IV. (0) based on the Minkowski space-time result, in the regime of field mass M ^> H, is 
corroborated by the derivation of a decay rate expression in de Sitter space-time shown in Ref. |32| , where it was shown 
that the decay rate behaves similar to the Minkowski one, but in a thermal bath at the Hawking temperature. However 
since in our calculation M x 3> H, this modification has negligible effect. There is also a simple way of understanding 
the result (|5.42|) in the case of an expanding space-time within the regime of parameters we are examining. Recall 
in conformal time (with conformal rescaled fields) our original model is no different from the one of a Minkowski 
space-time, except for the proper rescalings of dispersion relations and masses, u) = \/q 2 /a 2 + M 2 — > u> — \j q 2 + M 2 , 
where M = aM, Thus a rate evaluated in a conformally invariant theory, in conformal variables, is identical to that 
in flat space-time. For instance, consider the change of number of particles given in its simplest form, in conformal 
time, as 

dn - , 

— ~ Tn , (5.44) 

(XT 

which in terms of physical time becomes 

dn dr - f 

— ~ — rn = -n . (5.45) 

dt dt a y ' 

The above expression explicitly displays the (conformal) rate as suppressed by the scale factor. However take the 
case of, e.g., fermion production as given before, but in conformal rescaled quantities. Since the rate in conformal 
variables is identical to that of flat space-time, we have that 

■ ~ * ** h kj 8 - [ H ) ~ ^ h kj ^ { M i ) ' 

which then reproduces the above stated result Eq. (|5.42(l . 

We now substitute r^ = and T x given by (|5.42|) in l|5.41|l and consider the parameter regime relevant to our 
analysis here, where the couplings are A ~ O(10 -13 ) and g,g', h > O(10 -2 ). In this regime, we can generically drop 
the first nonlocal term in l|5.41|l . which comes from the $ scalar self-coupling, since this term is much smaller in 
magnitude than the dissipative term due to the \ corrections. Equivalcntly stated, since for the above parameter 
values the contribution to the dynamics of the background field ip coming from the <f> quantum modes are neglegible 
compared to those due to the Xj ones, we could as well consider from the beginning the original inflaton field $ 
in Eq. (|3.1|) as simply a classical (homogeneous) field <& = tp(t) in interaction with the remaining (quantum) fields 
in (|3.1|) . This was, for instance, the approach taken in Q- Adopting this approach, the only effect on the above 
calculations would be to drop all quantum inflaton self- interactions and thus keep only the nonlocal \ term in Eqs. 
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(|5.28|l and (|5.41|) . This point is important to note, since for the parameter regime to be considered here, we will 
see that while M x ^> H the same does not hold true for the inflaton mass, M^ ~ H, and so applying the WKB 
approximation for the <f> quantum modes could become questionable. However the smallness of the <fi self-coupling 
allows us to completely ignore the effects from the (f> quantum fluctuations. Thus, in what follows the quantum effects 
leading to the (^-effective EOM will only arise from the terms associated with the Xj dynamical quantum corrections 

in Eqs. (J5~2"£|) and JEHJ. 

Further approximations can be applied to (|5.41|) in the dynamical regime for which the motion of ip is slow. In 
this case an adiabatic-Markovian approximation can be applied to the nonlocal \- contributions. This converts (|5.41|) 
to one that is completely local in time, albeit with time derivative terms. The details of this approximation for 
Minkowski space-time can be found in Q|. Its extension to an expanding FRW background follows analogous lines. 
The Markovian approximation amounts to substituting t' — ► t in the arguments of the ^-fields in the second nonlocal 
term in Eq. (|5.41l) . The adiabatic approximation then requires self-consistently that all macroscopic motion is slow on 
the scale of microscopic motion, thus ip/(p, H < T x . Moreover when H < M x , the kernel K x (t, t 1 ) is well approximated 
by the nonexpanding limit H — * 0. The validity of all these approximations were examined in Q and they also will 
be examined in more detail in Sec. IVII The result of these approximations is that, after trivially integrating over the 
momentum integral in the last term in l|5.4I|l . the effective EOM Eq. (|5.41|) becomes [5J 

y+[3ff + T(^+ ^' fl) =0. (5.47) 

dtp 

By setting the couplings gj = g'- = g, h k j = h ~ g ^> A, the mass M x ~ gip >• m$ k and T x ~ N^h 2 M x /{STTUJ x ], it 
leads to the friction coefficient T(ip) in (|5.47|) 

TM = N X ^ a f (5.48) 

64ttM x Jl + a x <Ul + a x + 1 

where a x = N^h 2 /(8tt). 

In terms of the approximations used to derive (|5.47|) the equations (|5.35|) and H5.36[) also simplify. In particular, 
from (|5.3fcip) . we obtain that 

Pr ^ ?(<P)^j ■ (5.49) 

VI. NUMERICAL ANALYSIS OF EOM 

This section examines numerical results obtained from the basic equations that have been derived in the previous 
sections. In particular, we follow the considerations taken in the previous section as regarding the contribution of 
the nonlocal terms in the ^-effective EOM, dropping the neglegible 4> quantum corrections and keeping only the 
leading correction to the dynamics, given by the Xj nonlocal term. The behavior of the dissipative kernel coming 
from the Xj nonlocal term and solutions to the (^-effective EOM in the various approximations are then determined. 
Also, radiation production is studied in the different approximations. Finally the adiabatic conditions underlying the 
self-consistency of the basic equations in this paper are examined. 

A. The Dissipation Kernel 

Our analysis starts with the dissipative kernel. From D x (t, t') given in Eq. (|5.14l) and specializing our computations 
to the case of a de Sitter metric as appropriate for describing the inflationary phase, we can directly perform the time 
integrals appearing in D x (t, t') and then use the resulting expression in the dissipative kernel K x (t, t'), Eq. 15.20[) . to 
obtain 

K ,, ■/, f d 3 q u; x (c L a(t),t)co S [2W x (c L a(t),t,t')}+r x (c i a(t),t)sm[2W x (c i a(t),t,t')} 

K ^' > ~ J ( 27 r)3 w x (qo(t),tK(qa(*),«') [«5(qa(t),*) +I^(qa(t),t)] ^W°W'*'' > ^J 
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where 



W x (q,t u t 2 ) = 



dt"cJ x (q,t") = -jj<U! x (q,t)-U! x (t') + -^W: x 



2(6£-l)H 2 



x In 



u> x (t) 



2(6£-l)H 2 



u x (f) - Jm x + 2(6f - 1)H 2 




»'f 



2(6£-l)ff 2 



m 2 +2(6£-l).ff 2 



•2) 



and 



^x(<l'*i>*2) = exp 



-2m x r Xj (0) / dt"l/uj x (t") 



u> x (t) + Jm x + 2(6£ - 1)H 2 


u x (t>)-^m x + 2(6t;-l)H 2 


w x (t)-^m 2 x + 2^-l)H 2 


u x (t') + ^m x + 2(6Z-l)H 2 



-T x (0)/H 



(6.3) 



Note in arriving at Eq. H6.1|) . the factor of l/a(t) 3 in Eq. Ij5.14|l has been absorbed by a change of variable on the 
momentum integration. The solution for K x in Eq. (|6.1[) is valid up to errors of 0{u> /lo 2 ). It is useful to examine the 
behavior not only of K x (t, t') but also the integrated kernel 



I x (t)= / dt'K x (t',0). (6.4) 

Jo 

In Figs. ISJand El K x (t,0) and I x (t) respectively are plotted for the cases g — 0.1,0.37,0.4 and 0.5 in frames a-d 
respectively. In all the graphs the time interval 1/Ta(g) has been indicated, where To(g) = r x (0) as defined in Eq. 
(|5.43l) is the x decay width at zero momentum q = 0. The kernel K x (t,0) is seen to oscillate about zero with an 
overall enveloping amplitude that decays in a time interval ~ l/r x . The graphs of the integrated kernels I x {t), Fig. 
|HJ show that there is an overall skewness, and within the time interval of order 1/T X , the integrated kernels converge 
to almost constant values. Thus, although the kernel does not have a simple Gaussian or exponential decay behavior, 
the rapid oscillatory behavior that it does have effectively causes it to retain memory only over a time interval of order 
l/r x . It is also interesting to compare the kernel for the nonexpanding versus expanding cases, which is shown in 
Figs. 0and|S]for K x (t, 0) and I x (t) respectively at coupling g = 0.37. The graphs show very little difference between 
the two cases, which was expected since m x 3> H , and here is explicitly confirmed. For the parameters used in the 
figures we have that m x > 10 6 H . What differences there are between the nonexpanding and expanding space-time 
kernels become increasingly pronounced as t increases. This also is expected, since for very early times t <C l/H, the 
effect of expansion should be negligible. Note also in comparing Fig. ^Jwith Fig. [S]the y-axis is much more refined in 
the latter to help facilitate the desired comparison. However because of this in Fig. [S]the rapid decay of the integrated 
kernel below t < I/Tq can not be seen. 

From the graphs of I x (t) , the origin and validity of the local approximation Eq. I|5.47|l of the ^-evolution equation 
Eq. H5.28fl can be understood. For this, first note that irrespective of the effects that dissipative damping have on 
slowing the evolution of (p, a minimal damping always arises from the 3H<p term combined with the flatness of the 
potential, which in particular imply that within a time interval ~ 1/H, if and tp do not change significantly. In 
particular in integrating over the temporally nonlocal term in the yj-EOM Eq. I|5.28|l over a time interval of order 
1/H, ip and ip can be treated as constant and so taken out of the time integration. This leaves integration over only 
K x (t, 0) and as shown in Fig. El within a time interval ~ 1/T X , this integral rapidly converges to an almost constant 
value, which is precisely T/ip 2 of Eq. I|5.48|l . Since r x ^> H in all the cases in Fig. |S]it also means I x if) converges 
within a time t <C 1/H. 



B. (p effective equation of motion 



The solutions for ip(t) from the effective evolution equation are plotted in Fig. EJfor g = 0.1 and Fig. 1101 frames 
a-c for g — 0.37,0.4 and 0.5, respectively. For each case the solution is plotted from the exact one-loop evolution 
equation Eq. (|5.28|l (solid), the same nonlocal evolution equation except with the kernel being replaced with its 
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FIG. 5: The kernel K x (t, 0) of Eq. 16.11 for various interaction couplings g and h = g, with m x = gip(Q), <p(0) — JTipi, A = 10 13 , 
and £ = 0. 



nonexpanding space-time counterpart (dashed) , the adiabatic-Markovian evolution equation Eq. I|5.47|l (dot-dashed) 
and the evolution equation where no account for dissipative effects is treated (dotted). The latter dotted curves are 
the ones assumed in cold inflation studies, where the effects of dissipation are simply ignored. As Figs. El and HT)I 
indicate, this assumption can be critically wrong. In particular, for the large coupling cases in Fig. I1UI one sees that 
the effect of dissipation drastically affects the behavior of ip(t) from the underdamped evolution found in the dotted 
curves to overdamped evolution once dissipative effects are properly accounted for. Comparing the three curves in 
each frame which treat dissipative effects at different levels of approximation, we see that they are all in excellent 
agreement. In particular, the adiabatic-Markovian approximation, which is based on the simplified evolution equation 
Eq. (15.47(1 . is in excellent agreement with the exact evolution equation Eq. I|5.28|l . where the nonlocal kernel is fully 
treated numerically. Based on our examination of the kernels in Figs. ElandEl and the fact that the integrated kernels 
rapidly converge in a time t <C 1/H, these results for ip(t) come as no surprise. 

For the g = 0.37 and 0.4 cases, the longtime behavior appears to show oscillations as opposed to a complete 
overdamped relaxation. This is an artifact of our approximation of treating the ^-dependent \ field mass as fixed to 
the value of the field amplitude at the initial time to- This is done since then we can compute the kernel once and for 
all, before evolving ip in Eq. I|5.28|l . Allowing the \ mass to vary would require the kernel to be recomputed at every 
step of the evolution and that would be far too time consuming for this calculation be be tractable. However by doing 
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this simplification, it leads to the nonlocal damping term depending on the field amplitude as tp 2 whereas itshould 
be tp. This means as tp — > 0, our approximation causes the nonlocal term to go to zero faster than it actually should 
and in particular faster than the curvature of the potential, thus leading to the oscillations. Thus for g = 0.37 and 
0.4, the oscillations are simply an artifact of approximations used in numerically computing the tp effective evolution 
equation. 

Turning to Fig. Elfor g — 0.1, we find that the effect of dissipation is not significant enough to alter the evolution 
of tp(t) by very much. This is a weak dissipative regime, where in the adiabatic-Markovian approximation T < 3H. 
However the effect of dissipation is not entirely negligible. The two inset boxes in this figure close-up on tp(t) at the two 
extrema. They show that the amplitude of tp(t) is slightly less in the cases where dissipation is treated in comparison 
with the no dissipation (dotted curve) case. This indicates that energy is being depleted from the 92-system into 
radiation. 
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C. Radiation production 

In particular, Fig. 1111 shows the ratio /v /-ff for the cases g = 0.1,0.37,0.4 and 0.5 in a-d respectively. In each 
frame there is a plot of p r from the exact evolution equations, Eqs. I|5.28|l and (|5.36() (solid) and based on the 
adiabatic-Markovian evolution equations, Eqs. I|5.47|l and ()5.49(l (dashed). In all cases, the results show the exact and 
adiabatic-Markovian approximation are in good agreement. In particular, not only in the strong dissipative regime 
T > ZH but also in the weak dissipative regime T < 3iJ, the simple formula Eq. I|5.49|l for determining radiation 
production is valid. 
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D. Equilibration and thermalization: the asymptotic long time behavior 

So far we have not discussed the long time behavior of our results, in particular the equilibration and thermalization 
of the radiation produced by the dissipation mechanism discussed in this paper. Before doing that it is useful to show 
whether and how our results, particularly the numerical ones given above for both the dissipative kernel and for 
the evolution of the <p background field, can compare to recent numerical results obtained from the full evolution 
of fields far-from-equilibrium |33l | and therefore not restricted to quasi-equilibrium conditions only. Though here we 
have restricted our study of the dynamics for field configurations close to equilibrium and that evolve adiabatically, 
as we will see, our study still is able to capture many characteristics observed in the recent studies of the dynamics 
of scalar fields. 

Recently the authors in 33] have shown extensive numerical solutions for the kinetic (two-point correlation) equa- 
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FIG. 11: Evolution of radiation density p r , plotted as the ratio R = p r / H from various approximations and various interaction 
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tions for a scalar field in 1+1 D. Since these results seem also qualitatively to apply to 3+1 D, it is useful to see 
whether there is any similarity with the general characteristics observed for the dynamics, for both kernel and back- 
ground field, obtained here compared to those obtained in |33| . In particular, the authors of 33] have shown that 
the dynamics of correlations (that also applies to the time dependent number density evolution) can generically be 
divided into three basic regimes: a damping regime characterized by an exponential suppression of the correlations in 
a time scale idamp ~ 1/I\ followed by a drifting like evolution behavior characterized by smooth and slow changing 
of the modes, which typically lasts much longer than the initial damping evolution, after which a last regime sets 
in, the thermalization itself, in which thermal equilibrium is achieved, within a time scale ^thermal 3> id. 



amp- 



in the 

kinetic (or Boltzmann-like) approach in which the time evolution of correlations are solved, thermalization is seen as 
a direct consequence of self-consistently including scattering processes in the kinetic equations 33, 34, 35|. In our 
case, this would be equivalent to self-consistently take into account in our evolution equation and in the derivation of 
the nonequilibrium propagators, the backreaction of the produced radiation. This is fundamental in order to describe 
the thermalization process, since in this way proper equipartition of energy among the modes is taken into account 
and that will then lead to thermal equilibration in the long time evolution of the system field. Therefore, our results 
will not account for the very long time thermalization regime. On the other hand, we can check from our numerical 
results, in particular for the temporal behavior for both the nonlocal dissipative kernel and also for its time integrated 
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I|4.25|l . for various interaction couplings g = 0.1,0.37,0.4 and 0.5, with 



form, Figs. and El respectively, that they show similar behavior to the correlations obtained by the authors in Ref. 
[33j. In particular Figs. [5] and show for the kernel a quick damping of oscillations within the time scale 1/T X , 
which is then followed by a drift like behavior of almost unperturbed, small amplitude oscillatory like evolution over a 
much larger time scale than 1/T X . Up to the time scales we have studied the dynamics, we have seen no appreciable 
change in this behavior. Thus the way our kernel evolves with these two regimes is analogous to those observed in the 
full kinetic equations approach of [33| - This gives us an indication that, even though backreaction due to radiation 
production is not being fully taken into account, we are still capturing the relevant dynamics from the initial time of 
evolution (of no radiation) up to some very long time scale. Furthermore, it can be checked from the results shown 
in Fig. II II that during the time scale of evolution that we have studied, the produced radiation maintains a level that 
is only a fraction of the inflaton energy density, p r j ' p v < 10~ 5 , and, therefore, we expect its overall effect back on the 
evolution of the inflaton field to be only marginal. 

The full inclusion of scattering and then the description of the final, asymptotic equilibration and thermalization 
regime, could in principle be done within the kinetic, Boltzmann-like approach of 33, 34| l35| . or within other equivalent 
approaches able to describe the thermalization and equilibration process such as [3a, |37]]. The extension of these 
approaches to our expanding space-time multifield setting is an interesting (and likely most more complex) avenue 
for future work, but it is beyond the scope of this paper. 



E. Adiabatic and WKB approximations 



Finally we come to the analysis of some of our basic approximations used to derive the yj-EOM in its different forms. 
In Fig. ^]the validity of the adiabatic approximation is examined. From Eq. 14.25J I. recall this approximation is 
valid when H/u> + u>/u> 2 <C 1 (y-axis in the figure). For g = 0.5, Fig. HUIshows that ip remains overdamped throughout 
evolution and correspondingly Fig. 1 121 shows the adiabatic approximation remains excellent. However for g = 0.1, 0.37 
and 0.4 there are peaks crossing above one in Fig. 1121 thus meaning the adiabatic approximation breaks down in 
those regions. In comparing to Figs. ISland HUI all these peaks correspond to when the tp(t) evolution seizes to be 
overdamped and it goes through a maxima or minima. In these underdamped regimes, in any event the dissipative 
term has no significant influence on the evolution of the system and so the breakdown of this approximation is of little 
consequence. Moreover, in the context of inflation, for g = 0.37 and 0.4 these breakdown regimes first occur at very 
late stages near the end of inflation, and so are not in a regime of interest for large scale structure formation. Also as 
commented earlier, in order to make our calculation tractable on the computer, we treated the ^-dependent x field 
mass as fixed to the value of the field amplitude at the initial time to ■ This was done since then we can compute the 
kernel once and for all, before evolving <p in Eq. i|5.28[) . and thus cutting computation time by well over an order of 
magnitude and so bringing the computation time in the range of days as opposed to weeks. However by doing this 
simplification, it leads to the nonlocal damping term depending on the field amplitude as <p 2 whereas it should be ip. 
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This means as ip — ► 0, our approximation causes the nonlocal term to go to zero faster than it actually should and in 
particular faster than the curvature of the potential, thus leading to the oscillations. Thus for g = 0.37 and 0.4, the 
oscillations are simply an artifact of approximations used in numerically computing the <p effective evolution equation. 
For g = 0.1 the breakdown regimes of the adiabatic-Markovian approximation are real, however at early times, t < 50 
in the units shown in Fig. 1121 this approximation is excellent, so, for instance, results for radiation production at this 
time are reliable. 

VII. DISSIPATIVE MECHANISM IN SUPERSYMMETRY MODELS 

In the regimes where our dissipative mechanism is large enough to affect inflation, the interaction couplings also 
are significantly large to yield radiative corrections that harm the flatness of the inflaton effective potential. As such 
supersymmetry is needed, since it can cancel temporally local radiative effects from Bose and Fermi sectors, thus 
almost completely preserve the tree level potential. On the other hand, temporally nonlocal radiative effects, such as 
those that lead to dissipation, have very different space-time structure between Bose and Fermi sectors, and so are 
not canceled by SUSY. 

As discussed in Sec. IV CI the ^ x -fermions were included to mimic the effect of SUSY by cancelling the quantum 
corrections from the x _ bosons. However the basic dissipative mechanism we have been studying in this paper, light 
boson (inflaton) — » heavy boson — > light fermions, can be realized in very simple SUSY models. For example, in Ref. 
|5| we proposed the following model of two superfields 4> and X, 



W = -\/A$ 3 + g<PX 2 + 4mY 2 , 
o 



(7.1) 



where $ = <p + tp6 + 9 2 F and X = x + Sip x + 9 2 F X are chiral superfields. The field will be identified as the inflaton 
in this model with <j> = ip -\- a and {(p) = tp. This is the simplest SUSY model in which the inflaton has a monomial 
potential, in this case 



V (<p) = j<p\ 



(7.2) 



and which includes the standard reheating interaction term to an additional boson g 2 (p 2 x 2 ■ When ip ^ there is a 
nonzero vacuum energy and so SUSY is broken. This manifests in the splitting of masses between the x an d ip x SUSY 
partners with in particular 
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o (9 2 ~ o^5V 2 + Vlmgip + Am 2 
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= m l x - y^g^ 2 - 



(7.3) 



One can check that the one loop zero temperature effective potential correction in this case is not significant to alter 
the flatness of the tree level inflaton potential, 
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« V (<p) - -y. 



(7.4) 



The authors of Ref. |38| have recently studied independently the corrections in the model (|7.1|) including also the 
effect of finite temperature and reached analogous conclusion for the quantum corrections, that the T = and now 
also the thermal corrections can be kept under control. On the other hand this model has the interaction structure of 
the form Eq. 1)1 .2|) . and so leads to the dissipative mechanism studied in this paper. In particular, one of the Yukawa 
couplings of this model is Agxi^x^- Noting the mass splittings in Eqs. I|7.3|) . it means that the heavier x boson, xi 
can decay into a ip x fermion and an effectively massless inflatino ip. There will be a phase space suppression in this 
process due to the closeness in masses of ip x and xi so that the decay width now is 
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and this leads to the dissipative coefhcient being 

^ vVA 



2567r 2 m x 



V- (7.6) 



So in this case in general r xi ^^, x: ^, T < H. The radiation level, 1Z, during inflation is found from Eq. IJ5.49JI to be 

1/4 3/4 / \ 7/4 



^- V ^ °- 03 avH vJ • (7 - 7) 

Thus for A = 10 -13 and (p = mpi, 7?. > 1 arises for g > 0.73. However for 99 much larger than mpi, 1Z > 1 requires 
g > 1. Thus this model is not very robust in producing radiation during inflation, but nevertheless the effect also is 
not negligible. 

The radiation production during inflation in the model Eq. <|7.1|) can be greatly enhanced by adding some light 
fermions into which the Xi-bosons can decay. In any event, in a realistic particle physics model the inflaton sector, 
such as the model Eq. I|7.1|l . would interact with other fields. This could be done for example with an another 
superfield Y added to the superpotential Eq. (f T. l|t as hXYY/2, which leads to the Yukawa interaction term hx4>YipY- 
Provided Iray <C m x , the x decay width is unsuppressed and in particular would be just Eq. 15.42fl with all other 
subsequent expressions there also applicable here. For this model we find in the strong dissipative regime T > 3H 
that 



1/4 1 / \ 5/4 

ft S ^»2.66 j7 * , f^ • (7, 



Pr _ „ „„ 1 (mpi 

H ~ 5 3/4^1/2 A l/8 ^ ^ 

For A = 10~ 13 and if = 5mp\, which in cold inflation analysis of this model would be approximately where the 60th 
e-fold of inflation occurs, we get 1Z > 1 for g 3 / 2 h > 0.09. In the weak dissipative regime T < 3H we find 

1l= — ~ 0.036^^ [— j , (7.9) 

for which 1Z > 1 and the weak dissipation condition hold in the regime 0.00004 < g 3 / 2 h < 0.09. Thus warm inflation 
is very robust in this model. 

VIII. INFLUENCE OF DISSIPATION ON DENSITY PERTURBATIONS 

Provided that the radiation component present during inflation is bigger than the inflaton mass, p r > m^, one 
should generally expect that this radiation component will influence the fluctuations of the inflaton. Since the typical 

1 /4 

mass of the inflaton is ~ H , this amounts to the criteria already mentioned in the Introduction p r > H. Moreover, 
if one assumes thermalization, so that p r k, T, the inflaton fluctuations are then thermal. In this case the effect that 
the radiation component has on density fluctuations can be explicitly computed. Although it is beyond the scope 
of this paper to address the issue of thermalization, as a reasonable guideline thermalization is expected provided 
the decay width T x > H. In this section, some examples of density perturbations during the warm inflation regime 
will be presented and the differences will be compared to the comparable results that would be obtained under the 
assumption of cold inflation. 

For either ground state or thermal fluctuations of the inflaton, the density perturbations are obtained by the same 
expression 39], 

S B = ~S<p. (8.1) 

5 <f 
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In the warm inflation regime the fluctuations of the inflaton go in the strong dissipative regime as [I(j 

Sp 2 = (-\ y/HTT , for warm inflation (T > 3H), T > m , (8.2) 

and in the the weak dissipative regime as |41| 



, 3tt 



1/2 

Sp 2 = I — I HT , for warm inflation (T < 3H), T > m^. (8.3) 

In contrast, for cold inflation, where inflaton fluctuations are exclusively quantum |39| . 



H 2 

(2tt) 2 



Sp 2 — -{——To ; f° r cold inflation, T < m$. (8-4) 



The associated spectral indices, n s , for these three cases are |43 

d\nS 2 1/9 3 9 \ 
n s - 1 = " = - — e + -?y /3 , for warm inflation (T > 3H), T > m , (8.5) 



d\nk r \ 4 2 4 



17 3 1 
-— e + -ry — -/3 ) , for warm inflation (T < 3H), T > m^, (8.6) 



and 143 



n s — 1 = (— 6e + 2rj) , for cold inflation, T < m<p, (8-7) 

where k is the wavenumber of the inflaton mode and the slow-roll parameters are defined as e = m 2 V' 2 /(IdnrV 2 ), 
n = m%V"/(8irV), and (3 = ™2 T '^'/( 87rTT /). 

The spectral index will now be compared between the warm and cold inflation cases for the A(/> 4 /4 potential. For 
cold inflation the result are well known to be n s — 1 = — 3/N% |43|], where N e denotes the number of e- folds of inflation. 
So at iV e = 60, for example, n s — 1 = —1/20, with the model parameters A = 8 x 10 -14 and tp§o = 4.37mpi. These 
parameters correspond to the amplitude of the density perturbations of Sh ~ 10~ 5 . 

Turning to the warm inflation case, we now consider this potential coupled to additional fields in the manner of the 
Lagrangian Eq. (|3.1|l and account for the effects of radiation on density perturbations given by Eqs. I|8.2|l and (|8.3I) . In 
this case for the same model parameters and this same value y>60 of the field amplitude, these thermal effects increase 
the density perturbation normalization in the strong dissipative regime to 1 x 10 -3 . Thus the effect of dissipation and 
radiation production during inflation lead to a noticeable change in the behavior of the inflaton and its fluctuations. 
We now readjust the model parameters to properly normalize the density perturbations to the same value as before 
so that at N e w 60, Sh ~ 10~ 5 . Nevertheless the spectral index will still differ. In particular normalizing the density 
perturbations as before requires the parameters in the strong dissipative regime to now be A ~ 10~ 17 and the spectral 
index at 60-efolds becomes n s — 1 = —1.5/N e . So for N e = 60 this implies in the strong dissipative warm inflation 
regime the A(/(> 4 /4 potential leads to n s — 1 = 0.025 which is half the size of the correspond cold inflation result. Since 
the recent CMB satellite experiments, WMAP and upcoming Planck, should be able to discriminate spectral indices 
at the one percent level, the difference between warm and cold inflation might be detectable. This Section was simply 
illustrating some points regarding density perturbation differences in warm versus cold inflation. A detailed analysis 
of this issue will be presented in |44| . 

IX. CONCLUSION 

In this paper we have developed a formalism for treating dissipation in quantum field theory models with slowly 
evolving backgrounds in an expanding spacetime. The key steps for doing this were first computing the real-time 
matrix of dressed expanding space-time two-point Green's functions for the respective quantum fields in our system. 
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The solution of these Green's functions was obtained in a WKB approximation, which is valid for slow moving 
evolution of the background fields and for fields with masses much bigger than the Hubble scale. Having derived 
these Green's functions, we then used a standard response theory approximation approach for the derivation of the 
field averages appearing in the effective evolution equation for the background component of a scalar field <p = ($). 
The integration of the quantum field fluctuations employed a nonperturbative resummation. The resulting effective 
evolution equation for the background field <p showed dissipative features. 

As seen, our dissipative formalism differs from those used in treatments of reheating after inflation. In those 
cases, one is studying a fast moving, oscillating background component, typically in a linear relaxation (small field 
amplitude) regime. In contrast our analysis is applicable for slowly moving background fields that do not oscillate and 
in the nonlinear regime for the system field (the inflaton). As shown in Sec. |H1 the basic physics that underlies the 
dissipation in the two cases are markedly different. We have applied the nonlinear dissipative mechanism developed 
here to the inflationary regime. However, this same dissipation mechanism could also apply to preheating scenarios (if 
they are allowed by the model and given set of parameters), where the linearized, perturbative approximation for the 
inflaton breaks down and nonlinear, nonperturbative effects, such as the one studied in this paper, become important. 

In addition to deriving the basic equations for our dissipative formalism, a detailed numerical analysis was done. 
In particular, the key quantity that our formalism determines is temporally nonlocal terms that must be included in 
the ip background field effective evolution equation. In Sec. IVII these nonlocal kernels were numerically calculated 
and compared at various levels of approximations. For instance, the ip effective EOM was computed with the exact 
one-loop expression in Eq. (|5.28|) . A key question was the regime of validity of the simplified adiabatic-Markovian 
approximation Eq. IJ5.47I) to the exact equation. In Figs. 151 and ITUI these comparisons were made. In the regime 
where the WKB self consistency conditions Eq. I|4.25|) . together with the condition T Xj > H are satisfied, we found 
in Figs. 151 and ITUI that the evolution equation computed from the exact one-loop expression and from the adiabatic- 
Markovian approximation agree very well within the region of parameters we have concentrated our study. We also 
checked in Sec. IVII the radiation production from the ip system that emerges through dissipation. Once again the 
exact numerical treatment and the adiabatic-Markovian approximation agreed very well in the same regimes as for 
the evolution equation. These results are of great practical use, since calculating the exact numerical solution to 
the effective evolution equation is very time consuming on the computer, whereas the evolution equation in the 
adiabatic-Markovian approximation can be analyzed analytically. 

The immediate application of our dissipative formalism is to inflationary cosmology, in particular to determine 
warm inflation regimes and their properties. In Sees. IVIII and IVIIII we examined some consequences. In general there 
are two sorts of qualitative effects that dissipation can have on the inflationary phase. First, the evolution of the 
background field can be altered due to the nonlocal terms. The most dramatic example of that was shown in Fig. 
IIUI where accounting for dissipation, ip evolution was overdamped, whereas if one simply ignored these effects, the 
evolution would have the underdamped oscillatory behavior typically assumed in cold inflation studies. In particular, 
in the larger perturbative coupling regimes studied in Fig. IIUI the inflaton field would never have a reheating phase (in 
the sense of a fast oscillatory like regime for <p). It would simply relax to the minimum of the potential monotonically, 
dissipating radiation along the way, thus ending the inflation phase and initiating a radiation dominated phase. Less 
dramatic to this, but down to much lower interaction couplings, even though the background inflaton field evolution 
is not noticeably altered from dissipative effects, radiation is still being produced during inflation from conversion 
of vacuum energy. Although there are detailed questions about thermalization, that are beyond the scope of this 
paper to address, as a reasonable criteria when p r > H during inflation, one should expect that this radiation 
component will influence the inflaton fluctuations significantly from its zero-temperature zero-point level. To gain a 
better understanding of the extent that this radiation can influence the inflaton fluctuations, and hence the primordial 
seeds of density fluctuations, in Sec. IVIIII we computed the density perturbations when the radiation component is 
accounted for, and compared that to the naive expectation when the fluctuations are assumed to be zero-point ground 
state fluctuations. 

The key result of this study when applied to inflation has been that dissipative effects are predicted to occur during 
inflation in typical inflation models. These effects alter the single picture of inflationary dynamics assumed up to now, 
which we call cold inflation, into another possibility which we call warm inflation. To make accurate predictions from 
inflation models, which is now required for current high-precision CMB measurements, these dissipative effects must 
be treated. Moreover, dissipation effects can lead to some attractive theoretical consequences in inflation models. 
For example, for those parameter regions feasible to inflation and where the nonlinear and nonperturbative effects we 
studied here can become important, the emergence of effective strong dissipative phenomena are able, for instance, 
to sustain and drive inflation longer than when these dynamical effects are neglected [45j. Several other results also 
can follow as a consequence of the dissipative regimes studied here. In particular various problems, namely 77 |46J . 
graceful exit |j , quantum-to-classical transition [3, E3 , large inflaton amplitude |4(j , and aspects of initial conditions 
|48l l49| . can be remedied simply by properly accounting for the dissipative effects already in the model, rather than 
relying on additional modifications to the model, as is often done. 
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APPENDIX A: THE EFFECTIVE POTENTIAL AND RENORMALIZATION 

Consider the local terms (|5.25() appearing in the (^-effective EOM and associated to the field derivative of the 
effective potential for ip, 
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where we used Eqs. (|5.23|) and i|5.24f) . The momentum integrals in Eq. i|Alf) are divergent and require appropriate 
renormalization, which we perform here just as in standard Minkowski space-time. In Eq. ijAlf) we add mass and 
couplings renormalization counterterms to the classical potential so as 
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and we will consider from now on the masses and coupling constants as being the renormalized ones. The counterterms 
SrrA, <5£ and SX are fixed by the choice of renormalization conditions |2rll30j) 
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where we have chosen renormalization points ip — fx v and R — (jlj> in the above conditions so that the results are 
infrared finite in the limit of vanishing masses m^ and m x . Of course, these renormalization points are completely 
arbitrary and related to different choices by the corresponding renormalization group equations. 

Using an upper momentum cutoff A and changing to the physical momentum k p = k/a and cutoff A p = A/a, the 
momentum integrals in Eq. (jAlf) are easily evaluated leading to (for A — » 00) 
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We can now use the renormalization conditions (|A3|) in (|A4|) leading, for massless bare fields m^ = m x . = 0, to the 
rcnormalized expression 
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We can also extend the result for the renormalized effective potential when there is an additional coupling of $ to 
fcrmions ip x , in which case there is the additional contribution to IjAlfl 
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which leads to the additional contribution to l|A5|l (for bare massless ip x fermions, m^ = 0) 
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